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ABSTRACT 


A  representative  foimvilation  of  the  one  strip  hlunt  body  integral  method 
is  treated  in  svifficient  detail  to  delineate  the  numerical  difficulties 
involved  in  its  mechanization  as  well  as  to  identify  limitations  in  its 
application.  While  the  integration  of  the  governing  system  of  equations 
for  general  body  shapes  comprises  a  two-point  boundary  value  problem, 
the  partic\ilar  case  of  a  flat-faced  cylinder  is  shown  to  reduce  to  a  pro¬ 
cedure  involving  a  single  integration.  An  investigation  is  made  of  the 
extent  to  which  conservation  laws  are  satisfied  throughout  the  shock 
layer  region.  Significant  deviations  are  found  to  exist  depending  iqxjn 
the  method  used  for  computing  the  distribution  of  flow  properties  across 
the  shock  layer.  A  procedure  has  been  chosen  which  insures  compatible 
tangential  and  normal  flow  property  gradients  at  the  shock  and  body  while 
significantly  reducing  mass  and  momentm  defects  in  most  Instances .  Com¬ 
puted  flow  field  data  are  presented  for  spheres  and  flat-faced  cylinders. 
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Symbol 

a 

h 

i 

M 

P 

V 

Ar 

<r 

s 

£* 
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Speed  of  sound 
Enthalpy 

Mach  number 
Pressvire 

Body  radlxxs  of  curvature 
Shock  radixis  of  curvature 
Gas  constant 
Polar  coordinates 

Body  oriented  curvilinear  coordinates 
Entropy 

Velocity  component 
Total  velocity 
Cylindrical  coordinates 
Ratio  of  specific  heats, 

Distance  fron  hody  to  shock  measured  nonnal  to 
body  sxjrface 

Shock  detachment  distance  at  axis  of  symmetry 

Ancle  between  tangent  to  body  and  axis  of  symmetry 

Angle  between  tangent  to  streamline  and  axis  of 
symmetry 
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NOMBICLATURE 


Symbol 

Mach  angle, 
Density 


Angle  betveen  tcuigent  to  shock  and  vertical  axis 
Stream  function 


Subscripts 

sS 

n 

X 

o 

ao 


Denotes  ccmponent  In  <5  direction 
Denotes  ccmponent  in  direction 
Denotes  ccmponent  in  X  direction 
Quantity  evaluated  on  body  surface 
Quantity  evaluated  on  shock  wave 
Free  stream  valvie 
Stagnation  value 


Ihe  remaining  symbols  and  subscripts  are  defined  in  the  text. 
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1.  INTRODUCTION 


Die  Dorodnltsyn  integral  method  and  its  specific  applications  to  mixed 
flow  problems  are  well  represented  in  the  literature,  e.g.,  Refs.  1 
throu^  7.  Owing  to  the  nonlinear  natiire  of  the  two -point  boundary  value 
problem  involved,  however,  mechanization  of  the  method  involves  computing 
problems  of  a  nonroutine  nature.  The  present  report  discusses  these  prob¬ 
lems  in  detail  as  they  arise  in  the  (Belotserkovskii )  one  strip  approxi¬ 
mation.  Further  resxilts  concerning  the  basic  techniques  involved  as  well 
as  applications  to  blunt-nosed  configurations  are  presented.  Particular 
attention  is  paid  to  the  determination  of  flow  properties  across  the 
shock  layer.  For  the  sake  of  completeness  the  development  of  equations 
is  presented  in  detail.  Uie  procedure  and  attendant  notation  largely 
follow  Ref.  5. 


2.  DEVELOFMEMT  OF  EQUATIONS 


nie  equations  expressing  conservation  of  mass,  momentum  and  energy  for  a 
perfect,  Invlscld,  non -heat  conducting  gas  are,  respectively, 

(2.1) 

-o  .  (2.3) 

i*ere  ^  a  _ ^  ^ 

W  7i^  f 


Velocity,  pressure  and  density  have  been  nondlmenslonallzed  throu^  use 
of  maxlmxm  velocity  and  free  stream  stagnation  pressure  and  density, 
nils  choice  of  parameters  has  particular  advantage  >dien  applied  to  the 
Bernoulli  equation. 


For  general,  curvilinear  orthogonal  coordinates  where 

,  the 


vector  operations  are  given  by 
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Since  only  axlsynmetrlc  and  tvo -dimensional  configurations  will  be  treated, 
a  body-oriented  curvilinear  coordinate  system  (Fig.  l)  is  defined  such 
that 


FIGURE  1 


Substituting  the  Bernoulli  equation  for  the  mcmentm  equation  in  the 


J-direction,  the  governing  system  of  equations  hecanes 

(2-M 

«  O 

(2-5) 

\ 

II 

(2-6) 

(2-7) 

Let  ^‘Z-4 

(a) 

(D) 

/y-  jsiit-^^/o 

(c) 

(d) 

(2-8) 

(e) 

* 

'^‘^9  , 

(f) 

_L. 

,Aere 

Using  Bqs.  (2-6),  (2.7)  and  (2.9),  it  follows  that 

(2-10) 
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A  modified  continuity  equation  is  obtained  throu^  use  of  Eqs.  (2-8a), 
(2 -8b),  (2-10)  and  the  relation 


(2-11) 


In  addition,  Eqs.  (2-4)  and  (2-5)  are  ccmbined  to  obtain  a  modified  momen- 
tvm  eqxoation  in  "divergence"  form.* 


(2-12) 


Integrand  Ap-proacimations 

Consider  an  eqxiation  given  in  divergence  form,  i.e., 


*More  precisely,  one  takes  Eq.  (2 -5 ) /•  ‘^  **r  Eq .  (2-4)  and  adds  and 

subtracts 
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Integrating  along  an  arbitrary  line  in  the  /9-direction,  e.g.,  across  the 
shock  layer  frcm  the  body,  ,  to  the  shock,  ,  (llg.  2),  one 
obtains 


If  the  functions  are  approxijnated  in  the  form 

^  » 


it  follows  that 


4 
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Applying  Leibnitz's  njle,  one  obtains  an  ordinary,  first-order  differen¬ 
tial  equation: 


Ihe  integrand  approximations  employed  in  the  jo^sent  equations  are 


(2-13) 


Ihe  only  formal  mathematical  restriction  on  the  choice  of  representation 
is  that  the  ri^t  hand  side  of  Bqs.  (2-13)  te  a  specified  function  of  D 
containing  sf-dependent  parameters  vhlch  are  evaluated  at  either  the 
shock  or  body  (for  the  one  strip  approximation).  For  example,  one  may 
specify 


or 


In  some  instances,  the  specified  function  of  O  may  lead  to  behavior  in 
direct  contradiction  to  knovn  inviscid  results,  e.g.,  incorrect  surface 
or  shock  gradients  of  flow  properties.  A  more  complete  discussion  of  the 
validity  of  the  approximations  is  given  under  Inverse  Procedures. 


Ax^lying  the  preceding  considerations  to  the  modified  continuity  and 
momentvan  equations  yields  (noting  that  ^*.2^*^): 


(2-14) 


(2-15) 


1 
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Bqs.  (2-l4)  and  (2-15)  niay  be  expressed  entirely  in  terms  of  the  dependent 
variables  and  ^  throu^  use  of  the  following  auxiliary 

relationships : 


oblique 

shock 

relations 


'  T 

i 
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I 
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Expanding  Bq,s.  (2-l4)  and  (2-15)  and  introducing  a  geometrical  relationship 
between  tfT  and  ^ ,  one  obtains  the  governing  system  of  equations  in  final 
form: 


f 


(2-l6) 


(2-17) 


where 


(2-18) 


if  = 


r-..^^l(-y-4*&,-ry-^V%] 

yf.  -44^44'  ] 

f.  ^r'-4:f^i/-  ^4) 

4AiO-  a^Uy)  -^'Y\ 

/ _ z _ 

»  ^Orz-t/*) 
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3.  IMTEGRATION  PROCEDURE 
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In  order  to  initiate  the  numerical  integration  of  the  system  of  nonlinear 
total  differential  eqmtions,  knowledge  of  the  dependent  variables  ^ 
and  ^  at  the  axis  is  required.  Ihe  initial  values 

and  are  indeterminate  at  the  axis  and  require  special  treat¬ 
ment.  One  may  deduce  immediately  that  ;  however,  the 

shock  detachment  distance  ^  remains  unknown  euid  must  be  deter¬ 

mined  from  an  avKiliary  condition.  A  two-point  boundary  value  problem 
therefore  exists  which  must  be  solved  by  means  of  an  iterative  approach. 
(A  bar  will  be  \jsed  to  distinguish  the  "precise"  or  final  value  of  ^  as 
opposed  to  an  estimated  value  during  iteration.)  Bie  specific  approach 
taken  depends  on  nose  geometry  and  its  influence  on  the  surface  '^locity 
gradient . 


A  few  preliminary  remarks  are  in  order  before  discussing  the  two  types  of 
integration  procedures.  Application  of  the  Integral  method  to  smooth 
bodies  (i.e.,  having  continuous  slope)  is  facilitated  by  the  choice  of  a 
natijral  coordinate  system  such  that  the  body  surface  coincides  with  a 
coordinate  surface.  In  particular,  this  choice  insures  the  vanishing  of 
the  normal  component  of  the  surface  velocity  which  leads  to  certain  ana¬ 
lytical  simplifications.  Qnploying  body-oriented  orthogonal  curvilinear 
coordinates,  as  in  the  present  case,  one  is  only  practically  limited  by 
the  difficulty  in  analytically  describing  arc  length,  radius,  slope  and 
curvature  along  the  body  surface.  However,  when  a  slope  discontinuity 
(convex  comer)  is  encountered,  the  vj -derivatives  become  infinite  at  this 
point  and  an  "area  of  omission"  exists  which  prevents  continuation  of  the 
integration.  There  will  be  a  class  of  physically  meaningful  solvctions 
(for  the  region  from  the  axis  to  the  comer)  obtained  when  the  comer  cor¬ 
responds  to  the  body  sonic  point  (Plg.‘3)* 


II 


no  solvit  Ion 


0 

solvrtion  valid  for  \S  <vi" 

FIGURE  3 


A  mechanism  in  the  integral  method  \dvich  tends  to  further  limit  its  appli- 
cahility  in  certain  cases  is  the  direct  transmission  of  "body  Infonnation 
to  the  field  immediately  normal  to  the  body  rather  than,  for  example, 
along  a  characteristic  direction  in  the  supersonic  regime,  e.g..  Pig.  h. 


Zone  of  Influence 


FIGURE  U 
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Ohtis  a  change  in  body  slope  or  curvature  Imposes  an  immediate  perturbation 
across  the  adjoining  field.  This  effect  should  be  less  pronounced  in  the 
transonic  region  where  the  "normal"  mode  of  information  transmittal  is 
more  nearly  realized  in  the  physical  case.  It  may  be  reasoned  that  this 
effect  will  be  suppressed  as  the  integral  approximations  are  improved 
throu^  further  subdivision  of  the  shock  layer.  However,  by  rewriting  Eq. 
(2-l6)  (which  is  a  geometrical  relationship  independent  of  the  integral 
approximations)  as 


it  can  be  seen  that  a  curvature  discontinuity,  e.g.,  hemisphere-cylinder 
which  involves  a  change  in  curvat\ire  from  unity  to  zero,  results  in  a 
discontinuity  in  at  the  same  value  of  .  Further,  this  discon¬ 
tinuity  is  introduced  into  the  expressions  for  and 

(which  contain  ).  Caution  must  therefore  be  exercised  in  interpret¬ 

ing  flow  field  results  in  cases  involving  rapid  changes  in  body  geometry. 

3.1  Integration  Scheme;  Continuous  at 

Assuming  that  the  body  slope  is  defined  by  a  continiious  and  monotonically 
decreasing  function,  both  physical  and  mathematical  reasoning  calls  for 
continuity  of  the  surface  velocity  gradient.  The  denominator  of  the  ex¬ 
pression  for  f  however,  vanishes  identically  at  the  sonic  point. 

Hie  auxiliary  condition  for  determining  therefore  becomes 

in  order  that  finite.  Further  justification  for  this  proce- 

dure  is  given  in  Appendix  B. 

3.2  Integration  Scheme;  Discontinuous  at 

When  the  sonic  point  is  known  to  coincide  with  a  convex  comer,  a  dis¬ 
continuity  in  the  velocity  gradient  at  this  point  is  consistent  with 


13 


reqtilrements  of  inviscid  theory.  Hie  previovis  svDCiliary  condition  is 
therefore  no  longer  valid.  One  mvist  now  vary  €  tmtil 

^Ihere  is  a  class  of  "bodies  for  which  a  definitive  choice  of  method  cannot 
be  made,  e.g.,  sharp-cornered  sjJierical  caps.  For  example,  if  the  sonic 
point  occurs  below  a  convex  comer  (on  a  smooth  portion  of  the  body),  it 
will  not  be  possible  to  attain  the  requirement  J*  *  •  Further, 

if  the  corner  lies  below  the  intersection  of  the  "limiting  oh8n:^cteristic"* 
and  the  body,  the  slope  discontinuity  can  in  principle  still  affect  the 
upstream  subsonic  and  transonic  flow  regions.  In  practice,  however, 
there  is  no  feedback  scheme  in  the  one  strip  integral  method  to  properly 
accomit  for  this  effect. 

3.3  Initial  Values 


3*3*1  Axlsymmetric  Case 

Direct  substitution  of  the  initial  values  into  Eqs.  (2-l6),  (2-17)  and 
(2-l8)  resvilts  in  indeterminate  (o/o)  expressions  for  and 

^  .  Application  of  L'Hospital's  rvde  yields  the  proper  starting 

values. 


*Ref .  13  defines  the  limiting  characteristic  as  the  "locus  of  points  each 
of  which  has  only  one  point  of  the  sonic  line  in  its  zone  of  action." 
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vhere 


Also,  by  symnetry 


3.3-2  ^0 -Dimensional  Case 


No  indeteminatenesc  arises  in  this  case  (primarily  due  to  the  absence 
of  the  radius  terms);  therefore,  no  special  formulation  is  required.  It 
is  interesting  to  note  the  relationship  existing  betveen  two-dimensional 
and  axisymmetric  initial  derivatives,  i.e.. 


(3-1) 


(3-2) 
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4.  CONVERGENCE  SCHEME  -  SMOOTH  BODIES 


As  indicated  earlier,  the  determination  of  the  correct  initial  valvies  of 
^(TO)  involves  the  expression  for  the  surface  velocity  gradient  which 
for  convenience  may  be  expressed  as  /\/,  .  Assming  that 

N[^ i  f  L'Hospital's  rule  states  that  if  there 

ejcists  a  neighborhood  of  such  that  except  for  vf  , 

and  /V tiS*)  and  exist  and  do  not  vanish  simultaneously,  then 

whenever  the  limit  on  the  rl^t  exists.  How¬ 
ever,  in  the  present  case,  for  ^  ,  A/ {6 j  ^S*)  ^A/*^0  and 

j|C 

displays  singular  behavior  in  the  neighborhood  of  sf  .  One  must  also 
note  that  since  is  a  function  of  ,  it  is  in  effect  a  "floating" 
singiilar  point. 

Basically,  the  iterative  process  involves  successive  refinements  of 
such  that  /V^  ^  is  reduced  in  magnitude  to  the  extent  that  the  singular 
behavior  of  is  confined  to  a  small  nei^borhood  of  vT  .  At  this 

point,  extrapolation  techniques  may  be  applied  with  little  loss  in  accu- 
racy  in  order  to  continue  the  integration  beyond  sT^nT  . 

4.1  Initial  Value  Sensitivity  Stvidy 

In  order  to  establish  the  feasibility  of  a  given  prediction-correction 
scheme,  it  is  necessary  to  establish  the  sensitivity  of  the  integration 
to  the  choice  of  €  .  Frcm  the  standpoint  of  physical  measurements,  it 
is  reasonable  to  expect  that  a  carefully  conducted  experiment  wi]J.  yield 
shock  standoff  distances  (in  terms  of  nose  radii)  accurate  to  two  signifi¬ 
cant  figures .  In  the  present  integral  method,  for  a  given  A/f^ »  ^ #  inte¬ 
grand  approximation,  etc.,  there  exists  a  precise  value  of  ^ ^  ('vAilch, 
rounded  off  to  two  or  three  figures,  is  physically  acceptable)  which  will 
ca\ise  the  solution  to  be  sufficiently  well  bdiaved  in  the  sonic  region 
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such  that  the  integration  may  he  continued  into  the  supersonic  regime. 

The  degree  of  precision  called  for  was  found  to  be  in  excess  of  ei^t 
decimal  places.*  Figs.  5  throu^  7  explicitly  illustrate  the  behavior 
of  the  most  sensitive  function  involved,  i.e.,<J^^/^^  ,  aiid^O,  when 
£  is  perturbed  plus  or  minus  one  digit  frcm  a  reference  value  in  the 
second  throii^  eighth  decimal  places.  The  saddle  point  nature  of  the  sin¬ 
gularity  at  the  sonic  point  mentioned  by  previous  investigators  (e.g., 

Ref.  6)  is  clearly  evident.  Note  that  for  a  poor  guess  of  ^  ,  the  nei^- 
borhood  for  which  is  significantly  affected  extends  back  to  the 

stagnation  point,  e.g.,  for  ,  a  deviation  of  approximately 

six  percent  in  occvirs .  It  should  be  emphasized  that  the 

figures  presented  require  a  knowledge  of  ^  and  hence  represent  an 
"a  posteriori"  view  of  the  convergence  procedure. 

Initial  attempts  to  satisfy  the  downstream  sonic  condition  consisted  of 
a  series  of  "single  pass"  integrations  with  a  running  plot  being  kept  of 
A^^nd  'S*  in  order  t<^  make  successive  predictions  of  ^  (Pigs.  8  and 
9).  The  stopping  conditions  were  either  .  or 

Values  of  /V*"  and  *  were  determined  by  means  of  extrapolation  poly¬ 
nomials  (of  degree  applied  successively  to  ^  (to  determine 

vJ  )  and  /V"  [to  determine  A  positive  value  of  A/  indicated 

that  €  was  too  small  while  a  negative  value  indicated  that  €  was  too 
Isurge.  It  quickly  became  apparent  that  a  slope  prediction  technique 
alone  could  not  cope  with  the  nonlinear  bdiavior  of  the  functions  during 
early  stages  of  convergence. 

A  value  of  standoff  distance  established  for  A/  sufficiently  small  \rLll 
be  biased  by  choice  of  the  parameters  ^  and  ^ .  Therefore,  f  or  ^ 
are  to  be  changed  during  the  iterative  process,  a  high  degree  of  refine¬ 
ment  of  €  is  not  called  for  prior  to  the  change.  Regarding  the  paramEter, 


This  necessitated  the  luse  of  a  double-precision  mode  of  operation  (IBM 
7090-FORTRAN).  Uhfortunately,  double-precision  input  and  readout  capa¬ 
bilities  were  not  available;  therefore,  the  exact  nmber  of  figures 
reqvdred  beyond  el^t  could  not  be  established. 


€  SENSITIVITY  STUDY -N  VS  S 
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rji  one  must  obviously  compromise  between  fitting  from  too  far  (loss  of 
accuracy  in  extrapolation)  and  fitting  from  too  close  (loss  of  accuracy 
owing  to  erratic  behavior  of  functions).  Further  details  are  given  under 
Program  Degrees  of  Freedom. 


Ihe  convergence  procedure  employed  was  initiated  by  establishing  upper 
and  lower  limits  for  ^  (denoted  by  ^  and  respectively)  and  mak¬ 
ing  an  initial  guess  for  d"  *  ^  .  A  halving  mode  was  followed  until 

two  values  of  /y  were  recorded.  Fig.  10  displays  the  principal  opera¬ 
tions  involved.  Bie  tests  indicated  were  applied  after  each  integration 
step.  It  should  be  noted  that  the  curve  fitting  procedure  does  not  always 
yield  values  of  o  and  /V  ,  particularly  during  the  early  stages  of  con¬ 
vergence.  In  addition,  for  seme  cembinations  of  ^  and  step  size,  it  is 
not  possible  to  attain  a  value  of  \S  between  stations  corresponding  to 
^  bxA  ^ 

After  obtaining  two  values  of  /V ,  b.  linear  prediction  was  made : 


^  // 

If  >  d”  was  used;  otherwise  it  was  discarded  and  the  halv¬ 

ing  procedure  was  once  again  applied.  !Ihis  test  suppresses  the  "overshoot" 
tendency  of  the  slope  technique  in  the  nonlinear  range.  When  all  /V  h8id 
the  same  sign,  the  two  smallest  values  were  used.  If  both  positive  and 
negative  values  were  available,  the  smallest  value  of  each  was  used. 

Tirpically,  the  reduction  of  A/  proceeded  slowly  until  A/  ^  0{/0^) 
owing  to  the  fact  that  the  halving  technique  was  being  used  in  this  range. 
At  this  point  the  linear  prediction  came  into  play,  resioltlng  in  order  of 
magnitude  reductions  in  a  single  iteration,  e.g.,  decreases  in  A/  from 

and  frem  to  0^/0*^) for  a  given  itera¬ 

tion  have  been  noted.  No  particulsu:  advantage  appears  to  be  g^ned  from 


reducing  ,  \/  below  /O  .  Figure  11  displays  the  variation  of  the  ^ 
versus  'S  for 


No  attempt  was  made  to  vise  polynctnial  extrapolation  formulas  of  higher 
degree  for  estimating  new  values  of  ^ ,  One  disadvantage  of  increasing 
the  degree  is  the  associated  requirement  for  more  values  of  A/*,  i.e., 
polyncmial  of  degree  /r>  requires  /??*■/  determining  conditions. 


For  a  given  body  shape  and  specific  heat  ratio,  If  a  large  nvmber  of  runs 
corresponding  to  different  values  of  A/feo  "to  made,  one  may  expedite 
matters  throuj^i  use  of  polynomials  of  the  form 


<5^  « 


Z 


» 


using  the  first  few  cases  to  determine  the  constants.  Using  a  sixth 
degree  polyncmial,  it  >/as  found  that  five  reliable  figures  for  ^  could 
be  recovered,  therefore  reducing  the  number  of  iterations  required, 

4.2  Program  Degrees  of  R’eedcm 

Mechanization  of  the  integral  method  introduces  additional  parameters 
>Aiich  influence  the  numerical  resvats,  Hieir  values  are  set  primarily  by 
an  empirical  approach  using  the  compvtter  as  an  experimental  device.  A 
brief  discussion  of  their  influence  is  given  below.  The  data  are  taken 
frcm  sphere  calculations  • 

MM  ^ 

Integration  step  size;  Table  1  indicates  the  AS  effect  on  £"  and  v? 
for  a  fixed  stopping  condition  .SSCf*  •  The  fourth -order  Runge- 

Kutta  method  of  Integration  was  employed  in  ccmpiling  these  results . 

Use  of  other  numericeJ.  integration  schemes  was  not  considered  in  the 
present  investigation. 


25 


noun  II 


1 

I 

I 

I 

I 

r 

r 
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Meo  *  ^ 

.005 

.21592940 

.14823571 

.13692880 

.01 

.21592940 

.14823570 

.13692880 

.02 

.215929^^3 

.14823565 

.13692875 

.005 

.76711752 

.72000300 

.71176106 

s 

.01 

.76712370 

.71997798 

.71176337 

.02 

.76713051 

.71991991 

.71171965 

TABLE  1 

Near-sonic  stopping  condition;  For  those  cases  vhere  a  criterion 

Is  required  to  stop  the  Integration  procedure  short  of  the  sonic  point 
(defined  "by  JDC'S*)^0  ).  One  may  fix  this  point  In  terms  of  a  physical 
location  ^ S  or  may  specify  It  In  terms  of  a  certain  percentage,  ^  , 
of  the  sonic  velocity.  Althou^  the  former  method  requires  an  approximate 
knowledge  of  vJ*  ,  a  "linearizing"  effect  Is  Introduced  (see  Pigs.  12  and 
13)  with  reference  to  the  variation  ot  /V*  versus  .  Figs,  ik  throu^ 
l6  Illustrate  the  Influence  of  ^  on  the  standoff  distance  and  sonic 
point  location  as  well  as  the  quantities  ^  S  y 

c/J/^  evaluated  at  the  downstream  location  -J*  *  /  .  Ihe  latter  plots  are 
intended  to  indicate  the  increase  in  accuracy  gained  by  fitting  from  a 
point  closer  to  the  sonic  point.  An  asymptotic  bdaavior  is  Indicated  but 
computational  restrictions  prevent  determination  of  the  final  asymptotic 
value.  The  stopping  condition  ^  was  found  to  be  most  convenient 

from  a  compvrting  standpoint  in  spite  of  the  disadvantage  of  fitting  from 
different  terminal  values  of  vi*  [depending  iqxjn  the  value  of  standoff 
distance  chosen  which  in  turn  controls  the  functional  b^avior  of 
More  specifically,  a  randomness  is  introduced  in  the  fimction 
whidi  can  seriovisly  impair  slope  prediction  techniques.  In  addition,  a 
influence  is  introduced  unless  a  local  adjustment  of  step  size  is 
made  such  that  ^  as  precisely  at  the  end  of  a  given  step.  Table  2 
indicates  the  effect  of  on  ^  and 
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•Sf  variable  (DETERMINED  ON 
BASIS  OF  V|»  t  CRITERION) 


I 


nOUM  13 


7 

Mto*  6 

■ 

Co 

.975 

.21592946 

.14823571 

.13692880 

c 

.95 

.21592940 

.14823570 

.13692880 

.90 

.21592835 

.14823543 

.13692869 

.975 

.76717822 

.72006109 

.71182377 

.95 

.76712370 

.71997798 

.71176377 

.90 

.76682185 

.71952513 

.71135662 

SABLE  2 

Degree  of  curve  fit;  A  qtiadratic  expression  was  used  for  determining 
both  \S  and  A/  .  Other  investigators  (e.g.,  Ref.  7)  have  employed 
cubic  expressions.  A  parametric  study  of  the  influence  of  the  degree 
of  the  curve  was  not  made  for  smooth  bodies;  however,  results  pertaining 
to  this  topic  are  given  later  for  the  flat-faced  cylinder  case  (see 
Applications ) . 
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5.  COHVERGEKCE  SCHEME  -  BODIES  WITH  SONIC  CORKER 


As  stated  prevlovisly,  the  condition  to  he  satisfied  for  bodies  with  a 
single  convex  corner  at  which  sonic  velocity  is  reached  is  simply  given 
by 

Ihe  variation  of  with  €  for  a  typical  body  displays  none  of  the 
radically  nonlinear  behavior  of  the  smooth  body  case.  Qhrou^^  the  use 
of  a  quadratic  interpolation,  the  criterion  was 

satisfied  in  a  relatively  small  number  of  iterations.  For  the  special 
case  of  a  body  with  zero  cxirvatijre  up  to  the  comer,  i.e.,  flat-faced 
cylinder,  cone  with  detached  shock  wave,  no  iteration  is  required  in  a 
strict  sense  since  the  relation  is  precisely  linear.  Further  details 
are  given  under  Applications . 
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6.  INVERSE  PROCEDURE 


6 . 1  Iterative  Method 


Solution  of  the  governing  system  of  total  differential  equations  yields 
<!^)  )  their  derivatives.  Since 

and  have  "been  previously  given,  one  may  now  determine  ^  , 

functions  of  vj*  .  If  one  chooses 
to  follow  a  formal  inverse  procedure  in  calculating  shock  layer  proper¬ 
ties,  ^  ,  2^  ,  ^  ,  and  ^  must  be  evaluated.  Hence  , 

21^ySj/0)  and  are  specified  througji  use  of  Eq.  (2-13).  Finally 

/^ )  ^  }  ^  ^  determined  from  the  relations 


(6-1) 


€:i 


(6-2) 


(6-3) 


(6-4) 


Since  the  pressure  relationship  obtained  by  substituting  Bqs.  (6-2),  (6-3) 
and  (6-4)  into  Eq.  (6-l)  is  implicit  in  nature,  an  iterative  solution  is 
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required,  liifortunately,  the  solution  has  separate  roots  corresponding 
to  ^  i  and  switches  roots  along  normals  (  ^1'=  constant  lines)  cross¬ 
ing  the  sonic  line.  %e  calculated  flow  field  In  this  region  Is  always 
discontinuous  since  the  expression  does  not  allow  a  smooth  transition 
frcm  a  subsonic  body  point  to  a  supersonic  shock  point.  In  addition, 

for  the  case  of  a  flat -faced  cylinder  and  the  system  of  equations 
becomes  underspecified,  l.e.,  only  the  Z  and  linear  relations  plus 
Bernoulli  equation  are  available  for  determining  ^  .  Uils 

situation  also  arises  When  the  system  of  eqmtlons  Is  written  directly  In 
cylindrical  polar  coordinates  (e.g..  Ref.  8)  since  only  two  approximating 
functions  are  required. 

6.2  “linear**  Method 

Ref.  2  employs  separate  linear  approximations  for  and 

which  leave  the  differential  equations  unchanged  but  lead  to  a  simpler  but 

overspecified  algebraic  procedure  for  determining  the  flow  larpperties, 

l.e.. 


z’ fz’  ^  -^] 

4  - 


.Z 
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For  a  sphere  calculation,  this  situation  does  not  occur  if  polar  coordi¬ 
nates  3)  are  used  since  the  two  terms  collapse  natin^lly  to  a 
simpler  expression,  i.e,,  using  and  t  ^ 

reduces  to  .  Ref.  9  represents  S 

linearly  thxis  decoupling  the  p  tern.  Bi-s  decoupling  is  not  possible 
when  representing  linearly;  therefore,  while  the  original  gov¬ 
erning  equations  in  either  or  ^  coordinates  are  equivalent, 

the  integral  approximations  employed  are  not  equivalent. 

The  foregoing  considerations  indicate  that  use  of  a  formal  inverse  pro¬ 
cedure  in  the  present  one  strip  calculation  leads  to  ccmputational  dif- 
ficvilties.  Ref.  IT  (two  strip  approximation)  appears  to  circmvent  this 
difficulty  by  abandoning  the  formal  procedure  and  employing  qiadratic 
expressions  involving  values  of  the  flow  variables  on  the  body,  shock  and 
the  intermediate  center  strip. 

6.3  Gradient  Method 

derivatives  of  the  flow  variables  in  the  vf  and  /?  directions  at  the 
shock  ftwd  body  may  be  expressed  entirely  in  terms  of  quantities  yielded 
by  the  basic  integration  for  3 ,  ^  and  .  This  suggests  a  flow 
field  determination  based  vqxjn  a  knowledge  of  these  gradients  involving 
cubic  approximations  for  three  variables  with  the  remaining  variable 
being  recovered  from  the  Bernoulli  equation  and  the  entrojor  function. 

An  empiriceO.  procedure  involving  different  combinations  of  cubic  fits  led 
to  the  use  of  the  following  system  of  equations  (see  Appendix  I  for  the 
derivation  of  the  gradient  functions); 


*It  should  be  noted  that  an  approximate  higier  order  scheme  involving 
gradients  is  discussed  in  Ref.  5  in  idiidi  the  one  strip  differential 
equations  are  modified  with  correction  terms  In  an  attempt  to  obtain 
ijiq)roved  body  and  shock  data. 
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/)  f  ^  and  ^  are  determined  iising  expressions  of  the  form 


is  s\ibseq,uently  determined  from  the  Bernoulli  equation  with  the 
cubic  representation  (denoted  by  a  prime)  being  employed  in  order  to 
fix  the  sign. 


Also 


nhe  cdm  of  the  foregoing  procedxire  is  to  recover  maximun  Infoxmatlon  from 
those  elements  of  the  solution  \diidi  most  closely  agree  with  experimental 
data.  The  usefulness  of  the  gradient  method  will,  of  course,  be  dimin¬ 
ished  for  those  cases  for  Uhich  the  Aock  and  body  values  become  less 
accurate. 


39 


7.  CORSERVATION  REIATI0M5 


A  measure  of  the  mathematical  cozisistency  of  the  gradient  approach  is 
afforded  hy  detezmining  the  extent  to  ^diich  conservation  lavs  are  satis¬ 
fied  throuj^oxzt  the  shock  layer  region.  It  has  been  noted  (Ref.  10 )  that 
significant  deviations  in  conservation  of  mass  and  mcmentum  can  occur  in 
flov  fields  obtained  using  the  integral  method. 

Considering  the  control  voltxne  shown  in  Fig.  17,  the  conservation  equa¬ 
tions  are  given  by: 

Mass: 

/a) 

-  (^7r)Jjo4/^'ah  (7.1) 


Momentm: 


o 


(7-2) 


o 


(7-3) 


n 


For  an  adiabatic  flow  process  involving  a  thermally  €md  calorically  per¬ 
fect  gas,  {/}  nondijnensionalized  using  therefore, 

vising  Bq..  (2-7),  the  integral  fonn  of  the  energy  relation  reduces  to  the 
fom  of  the  mass  conservation  equation. 


(7-4) 


-  /  .  (7-5) 


values  of  HO  and  A/f  have  heen  ccmputed  for  a  ntanber  of  cases.  Die 
results  are  presented  In  Figs.  l8  throu^  20.  Ihe  deviations  noted  are 
primarily  a  conseqioence  of  the  inverse  procedure  employed  to  recover  the 
flow  variables. 


7.1  S-phere  Restilts 

Fig.  id  presents  a  conxiarlson  of  the  deviations  fi'ctn  mass  conservation 
resulting  frcm  application  of  the  "linear,”  iterative  and  gradient  inverse 
methods  for  A/f^=  3>  6,  9*  A  single  (  6)  curve  is  presented  for  the 

iterative  method.  Bie  break  in  this  curve  corresponds  to  the  "double¬ 
valued"  region  mentioned  previously.  Frcm  a  practical  standpoint,  the 
deviations  occurlng  in  the  range  corresponding  to  low  supersonic  Mach 
nvBibers  in  the  field  (.8  <sX  <’.9)  are  particularly  significant  in  that 
starting  solutions  for  characteristics  are  derived  from  the  flow  field 
results  in  this  region.  Stcppression  of  errors  which  could  subsequently 
propagate  into  the  downstream  field  is  therefore  necessary  in  this  range. 
In  this  respect,  the  gradient  method  provides  a  siibstantlal  is^rovement 
over  the  other  techniques.  Conservaticm  of  momentum  deviations  are  also 
significantly  reduced  through  use  of  the  gradient  method  (Fig.  19)* 

7.2  FLat-Ihced  Cylinder  Results 

As  indicated  previously,  the  need  for  two  of  the  approximating  functions 
vanishes  in  the  case  of  the  flat-faced  cylinder.  A  decision  was  made  to 
retain  these  relations  in  detemlnlng  the  flow  field  since  the  approxima¬ 
tions  consistent  with  this  "linear"  inverse  method  appear  to  have  more 
validity  than  the  gradient  method  on  the  basis  of  results  shown  in  Fig. 
20.  Since  the  accuracy  of  the  gradient  method  is  directly  related  to  the 
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- LINEAR 
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validity  of  shock  and  body  data,  particuleu*ly  "second-order”  quantities 
such  as  and  ,  the  deviations  in  interpreted  as 

an  indication  of  shortccoiings  in  the  basic  solution  for  this  body  shape 
and  at  the  given  free  stream  conditions. 
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8.  APPLICATIONS 


The  success  of  the  one  strip  Integral  method  in  predicting  shock  shapes 
and  surface  pressiire  dlstrlhutlons  has  heen  established  by  prevlovts  in¬ 
vestigators,  excepting  those  cases  \jhere  is  too  close  to  unity  or 
vhere  rapid  changes  in  body  curvature  are  Involved.  Ihfortunately,  little 
Infozmatlon  is  available  for  detemlnlng  the  accuracy  of  the  distribution 
of  flow  properties  across  the  shock  layer  [sonic  line  date  (Ref.  U)  and 
interferometric  ballistic  remge  data  (Ref.  12)]. 


In  view  of  the  preceding  considerations,  the  results  presented  do  not 
Include  experimental  data  but  are  intended  to  typify  the  information 
yielded  by  the  Integral  method. 


8.1  Sphere 

Ihe  variation  of  shock  detadiment  distance  with  free  stream  Mach  nvmber 
for  ^  -  1.4  is  given  in  Fig.  21.  A  sixth  degree  polyncmiea  of  the  form 


c_ 


was  fitted  to  the  nvinerlcal  data  in  order  to  provide  accxirate  detachment 
distances  for  intemediate  Mach  nmbers  (>dilch  are  usable  in  the  present 
formulation  to  approximately  five  decimal  places  as  Initial  €  values). 

A  few  cases  were  r\jn  for  for  1.2  and  A^«=  20, 

^  ^  .073666)  whldr  displayed  trends  in  accord  with  the  results  of  other 
theoretical  methods.  A  plot  of  sonic  point  location  versus  /V^  is  also 
inclirded  (Fig.  22).  For  4,  a  crossplot  indicated  that  the  relation- 

ship  between  ^  and  \S  beccmes  nearly  linear.  Fig.  23  presents  the 
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variation  of  shock  radius  of  curvature  at  the  axis  with  A/L .  Bae  limit- 

Om 

ing  value  of  as  is  significantly  higher  than  the  correspond¬ 

ing  result  obtained  by  assuming  a  concentric  shock  wave,  i.e.,/^ . 


The  distribvrtion  of  flow  properties  across  the  shock  layer  was  determined 
by  means  of  the  gradient  method.  In  particular,  the  variations  of 

and  are  shown  in  Pigs.  2k,  25  and  26,  respectively, 

for  6  and  1.4.  Die  density  plot  displays  the  small  jD  varia¬ 

tion  in  the  stagnation  region  which  is  the  physiceLL  basis  for  a  nunber  of 
approximate  constant  density  solutions. 


Pig.  27  gives  a  caapeurison  of  shock  shapes  derived  throu^  use  of  the 
method  of  characteristics  (with  an  integral  method  starting  line)  and 
the  integral  method  alone,  nie  characteristic  sepcurating  the  sphere 
(o  <  region  of  influence  and  the  region  additionally  influenced 

by  the  afterbody  is  indicated  in  order  to  illustrate  the  "illegal"  dis¬ 
turbance  imposed  in  the  integral  method  by  change  in  afterbody  gecmetry. 


Typically,  as  the  free  stream  Mach  nmber  is  increased,  the  integral 
method  shock  wave  tends  to  be  farther  below  the  true  shock  location  (as 
defined  by  characteristics  and  experiment)  for  -^>2  .  It  is  interest¬ 
ing  to  note  that  for  a  ccmplete  sphere  the  integration  into  the  super¬ 
sonic  regime  continues  smoothly  until  approaches  , 

point  the  nonnal  to  the  sphere  cannot  inter¬ 
sect  the  shock  wave.  For  example,  for  1000,  the  integration  pro¬ 

ceeded  up  to  \S  =  3«13;  i.e.,  179‘3®  around  the  sphere! 

8.2  Plat-P^ed  Cylinder 

!Ihe  parameters  which  control  a  given  integration  are 
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COMPARISON  OF  SHOCK  SHAPES  DERIVED  FROM  CHARACTERISTICS  AND  INTEGRAL  METHODS 


For  conical  or  wedge  nose  caps  teznlnatlng  In  a  sonic  comer 

the  latter  two  parameters  hecome  fixed  and 
hecane  directly  proportional  to  the  Independent  variable  As  a 

result,  no  coupling  exists  between  ^  and  the  geometrical  parameters  as 
the  integration  progresses.  Sie  location  of  the  body  sonic  point  which 
serves  as  the  only  characteristic  dimension  remaining  in  the  problem 
therefore  becomes  purely  a  function  of  S'  •  Hence,  a  scaling  is  Inq^led 
which  allows  one  to  fix  the  sonic  comer  radius,  ^  ,  at  a  pradiosen 
value  employing  a  maxlmun  of  two  Integrations.  Die  first  integration 
establishes  the  ratio 


and  the  second  integration  serves  to  adjust  S  *^0  a  convenient  ref- 
erence  value  using  ^  /fV  present  ease,  ^  was  chosen  as 

unity),  nie  latter  operation  was  found  to  be  more  convenient  than  one 
involving  a  storage  of  all  q,uantltles  to  be  sealed. 

Since  the  accuracy  of  the  Integral  method  is  directly  related  to  the 
approximations  governing  functional  babavlor  in  the  -direction,  the 
validity  of  employing  results  derived  from  local  transonic  solutions 
(Hef .  l6)  within  the  framework  of  the  computing  procedures  was  cheeked 
in  the  following  manner.  Ibe  sealing  property  of  the  flat-faced  cylinder 
ease  provided  a  means  for  detexnlnlng  the  consistency  of  using  the  func¬ 
tional  fom  suggested  by  second-order  transonic  theory  for  the  surface 
velocity  extrapolation  fomula  at  the  comer,  l.e.. 


« - 

A  refomulatlon  in  texns  of  the  coordinate  system  employed  is  necessaiT' 
to  handle  oases  for  Zn  this  instance,  sf  would  be  replaced 
by  a  corresponding  variable. 
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For  O  <r/7  <r  /  i  this  expression  eillows  the  upstream  surface  velocity 
gradient  to  approach  an  infinite  value  at  the  corner  (as  dictated  by 
inviscid  flov  theory). 


Hie  success  of  the  scaling  procedure  was  influenced  by  the  choice  of  the 
exponent  /7  .  For  example,  a  value  of  yielded  a  sonic  comer 

radius  of  1.0000001  while  a  value  of  /?  ■»  (given  by  second  order 
theory)  gave  a  "less  precise”  scaling,  1.0000044.  A  weak  depend¬ 

ence  upon  O  therefore  exists.  In  addition,  mmerical  detennination  of 
in  an  inverse  fashion  employing  the  last  valid  points  ccnputed  near  the 
comer  yielded  a  value  of  O.4996.  Finally,  employing  a  cubic  fit  which 
is  in  direct  contradiction  to  the  requirement  that  ^ 

as  O  gives  the  poorest  scaling,  -  1.0004.  S  C 


From  a  practical  standpoint,  the  preceding  argunents  become  somewhat 
academic  if  the  integration  is  not  continued  past  the  sonic  comer  since 
upstream  flow  field  results  are  not  significantly  affected  by  the  fitting 
proced\ires . 

Owing  to  the  steepening  of  the  velocity  gradient  near  the  comer,  a  step 
size  reduction  was  introduced  as  the  comer  was  approached.  Ihis  effect 
became  partlculsurly  critical  eis  the  Meudi  nvmber  was  decreased.  Ihe  reduc¬ 
tion  criterion  was  based  on  a  test  of  rather  than  ^  owing  to 

the  rapid  divergence  of  the  velocity  gradient  expression  within  the  spem 
of  a  few  integration  steps .  Sequential  interval  reductions  from  =  .01 
to^*S  =  .0001  (based  on  =  l.O)  were  employed  in  sane  cases. 

Qhe  scaling  relationship  which  exists  for  flat- faced  slabs  or  cylinders 
serves  an  auxiliary  purpose,  namely,  that  of  program  checkout  (excepting 
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those  terms  which  vanish  when  If  the  "two  integration  pro¬ 

cedure"  does  not  place  the  sonic  point  precisely  at  a  chosen  location, 
then  an  error  exists  in  the  program. 

8.3  Nose  Caps 

Having  determined  the  flow  field  for  a  smooth  nosed  configuration,  one 
may  proceed  to  generate  resvilts  for  a  family  of  nose  "caps"  provided  that 
they  teminate  in  a  sonic  comer.  Modification  of  extrapolation  proce¬ 
dures  near  the  sonic  point  is  necessary,  however,  to  allow  for  the  proper 
functional  behavior  in  this  nei($iborhood,  e.g.,  infinite  upstream 
gradient . 

Ihe  sphere  provides  a  useful  example  for  illustrating  the  preceding  con¬ 
siderations.  Representing  the  sphere  chock  standoff  distance  at  a  given 
and  ^  by  s',  each  integration  for  a  given  ^  will  produce 
results  corresponding  to  a  spherical  cap.  A  few  iterations  are  necessary 
in  order  to  determine  the  value  of  o  corresponding  to  a  specified  comer 
location)  however  the  convergence  is  relatively  rapid,  A  comparison  of 
results  for  complete  and  truncated  spheres  is  shown  in  Jig.  28.  Die 
Newtonian  value  is  also  included  for  comparative  purposes. 

8.U  Two-Dimensional  Results 

Plat -Paced  Plate 


Ptaploying  the  step  size  reduction  procedure  used  in  the  axisymraetric  case 
near  the  comer,  no  difficulty  was  encountered  in  verifying  the  C  -scaling 
property  for  the  two-dimensional  flat-faced  body,  i.e.,  only  two  integra¬ 
tions  were  required  to  fix  the  sonic  point  precisely  at  the  comer  of  the 
body.  Pigs.  29  and  30  give  a  comparison  of  shock  shapes  and  surface  pres¬ 
sure  and  velocity  distributions  for  the  flat -faced  plate  and  cylinder. 
Referring  to  Eq.  (3-2)  for 


se 


For  ,  therefore,  the  shock  curvature  at  the  axis  is  greater 

for  the  axisymmetric  case.  Cxis  inequality  was  satisfied  for  the  cases 

computed,  e.g.,  for  5*8>  2{.hg662)  ^.87025. 

€0 

Circular  Cylinder 

On  the  basis  of  preliminary  evidence,  the  convergence  procediare  initially 
devised  for  the  axisymmetric  case  was  inadequate  for  handling  two- 
dimensional  bodies,  i.e.,  circular  cylinder. 

Bxe  principal  factor  involved  was  the  singular  bdxavior  of  for 

the  same  range  of  ^  in  which  the  velocity  gradient  displays  erratic 
bdiavior.  Setting  tests  which  restricted  shock  curvature*  (see  Pig.  lO) 
resulted  in  successful  runs  at  discrete  Mach  nunbers  (e.g.,  5)» 

However,  the  iteration  displayed  generally  erratic  and  inconclusive  be¬ 
havior.  For  the  converged  runs,  the  shock  wave  displayed  points  of 
inflection  in  the  supersonic  regime  which  could  not  be  ascribed  to  known 
physical  behavior. 

It  is  worth  noting  that  the  USSR  cylinder  results  (Ref.  2)  were  presented 
primarily  for  the  two  and  three  strip  calculations  with  only  surface  pres¬ 
sure  data  given  for  the  one  strip  version  (indicating  the  possibility  of 
inadequacy  of  the  one  strip  results).  In  addition,  no  other  reference 
was  found  which  presented  two-dimensional  integral  method  results  of  any 
description. 


*For  the  axisymmetric  case.  Ref.  5  proposed  a  test  to  handle 

situations  where  the  initial  ^  guess  is  extremely  poor. 


9.  COMCLUDING  REMARKS 


In  the  present  report  emphasis  has  been  placed  on  the  discussion  of  a 
representative  formulation  of  the  "blunt  body  integral  method  which  in¬ 
cludes  sufficient  detail  to  delineate  the  numerical  pitfalls  involved 
as  well  as  some  practical  limitations  in  its  application.  Particular 
importance  was  placed  on  the  flow  field  determination  in  order  to  provide 
an  accvurate  description  of  the  flow  properties  across  the  shock  layer. 

Ihe  gradient  method  was  employed  (for  smooth  bodies)  since  its  applica¬ 
tion  Insured  compatible  tangential  and  normal  flow  property  gradients  at 
the  shock  and  body  while  minimizing  deviations  in  conservation  of  mass 
and  momentum  in  the  shock  layer. 

In  the  course  of  the  investigation,  a  previously  unreported  feature  of 
the  method  came  to  lle^t,  namely,  that  for  the  flat -faced  cylinder  case, 
the  two  i)oint  boundary  value  problem  reduces  to  one  involving  a  single 
integration  of  the  system  of  governing  equations.  Oils  featiire  obviously 
affords  a  significant  reduction  in  ccmputlng  time. 

Regarding  the  practical  utilization  of  the  method,  rational  flow  field 
data  have  been  generated  for  a  range  of  blunt-nosed  axisymmetrlc  configu¬ 
rations  employing  a  combination  of  the  Integral  method  and  the  method  of 
characteristics.  Application  of  the  latter  method  was  found  necessary 
in  order  to  maintain  downstream  accuracy. 
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APPENDIX  A 


FLOW  VARIABLE  GRADIENTS  AT  SHOCK  AND  BODY 


Ihe  governing  equations  expressed  in  body  oriented  curvilinear  coordinates 
are  given  by 


Shock  Gradients 


Uie  oblique  shock  eqxjations  serve  as  auxiliary  relations  in  obtaining  flow 
variable  gradients  immediately  behind  the  shock  wave .  In  addition,  geo¬ 
metrical  considerations  (Fig.  Al)  enable  one  to  express  -derivatives  in 
terms  of  /}-  and  ^-derivatives,  i.e.,  letting  ^ ^ ^ * 


fdA 


=  - _ 


'3.)  ^ 


or 


(A5) 
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Use  of  this  operator  in  Eqs.  (Al),  (A2),  (A3)  and  (A^)  results  In  an 
algebraic  syatem  of  four  equations  for  the  gradients  > 

stabsequently  obtains: 
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(A8) 


4(? 


#('^i)  ■ 


where 


-5^ 

y:^_  -/  /iA  <W  .  ^  <y'/^\  y 


a-  -/  2>  ^  4r^  ,  =  -J^;tkl 

^  .sw^  ^  ^(z-^^)  ^casfO 


(A9) 


4-- 


</V-/  » 


-  v: 


dA 


-  (cx^a^sA  *^v/3<^^ 


Ibe  data  required  to  evaluate  these  gradients  are 
‘%,^,J'  (note  that 


One  may  now  express  the  -derivatives  in  terms  of  the  /? -derivatives 
throuc^  use  of  the  original  equations. 


.IN 


(All) 

4*4 


’'4 

(A13) 
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[Alternative  expressions  may  be  determined  throvi£^  use  of  the  equations 
^  and 


Body  Gradients 

Ebcpressions  for  and  are  obtained  by  setting  0^0 

and  ^  M  O  in  the  governing  eqxaatlons. 


\diere  is  assvned  known.  Note  that  the  coefficient  of  the  sur¬ 

face  velocity  gradient  vanishes  at  the  sonic  point  indicating  that  similar 
behavior  in  the  integral  method  is  Independent  of  the  formulation  and 
apprcKlmatlons  enjoyed. 
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Die  remaining  normal  derivatives  are  determined  employing  a  procedure 
outlined  in  Ref.  15.  Consider  a  function  vith  and 

defined  as  coordinates  tangential  and  normal,  respectively,  to 
the  streamline  direction.  One  may  therefore  write 

^  ^  ^  ^ 

Assuming  that  (^=  shock  wave  angle  at  intersection  of 

shock  and  streamline  Z  -  constant),  the  expression  reduces  to 

2r 

SZ 


J 

1 

ir 


FIGURE  A2 

From  continuity  considerations 

-/ 

or 


^4  ^ 

r^' 
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Therefore 


Unally 


Ihe  tangential  derivatives  at  the  surface  are 


Letting  J*  represent  any  flow  variable,  the  angle  between  the  curve 
"f  =  constant  and  the  >i^-axis  is  given  by 


Hence,  this  angle  may  be  detemlned  at  the  shock  and  the  body  employing 
the  derivatives  given  above. 


APPENDK  B 


BEHAVIOR  OF  IHE  SURFACE  VELOCITY  GRADIEMT  EXPRESSION  IN  THE  MEIGHBORHOOD 
OF  THE  SONIC  POINT 


A  computing  scheme  was  devised  which  was  intended  to  eliminate  the  need 
for  curve  fitting  and  extensive  ^  refinement  required  to  cross  the  sonic 
point.  It  was  based  upon  the  following  assumptions: 

a.  Application  of  L' Hospital's  rule  to  the  equation  for  the  sxarface 
velocity  gradient  yields  a  well  behaved  expression,  exact  at 

•sS^'S*  for 

b.  Ihis  expression  approximates  the  correct  value  of  in 

some  finite  nei^borhood  of  'S  overlapping  the  region  in  which 
the  value  calculated  by  the  conventional  expression  is  still  valid. 

c.  Ihe  expression  is  less  sensitive  to  the  initial  value, ^  ,  than 
the  original  gradient  equation. 

Die  scheme  involved  switching  to  the  L'Hospital  or  "L"  value  Jtjst  before 

begins  to  behave  erratically  and  switching  back  on  the  other  side 
of  the  sonic  point,  e.g..  Fig.  Bl. 


PIGURE  Bl 
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Switching  expressions  resulted  in  a  temporary  instability  in  the  inte¬ 
gration  which  recovered  in  the  supersonic  region.  Characteristically, 

the  integration  process  displayed  strong  "recovery"  ability  in  that  an 

# 

erratic  value  introduced  by  stepping  too  close  to  si*  introduced  only 
a  temporary  perturbation  which  was  rapidly  damped  such  that  downstream 
results  were  relatively  unaffected. 

Mathematical  Considerations 

!nie  denominator  of  the  surface  velocity  gradient  equation  Is  given  by 


Khtr. 

for 

Now 

Therefore 

Assvnlng  that  the  surface  velocity  gradient  has  a  finite,  nonzero  value 
at  the  preceding  considerations  Indicate  that  application  of 
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L'Hospltal's  rule  yields  an  expression  \dilch  Is  well  'behaved  In  the  nel^- 
"borhood  of  the  sonic  point.  More  explicitly, 


///r? 

a(S 


o(S 


c/s  c/S 

cfS  ^ 

^  p 


(Bi^) 


d  -  ^  ^  d£  £,  y.  £- 

^  ^  9S  a<S  &<f  e/S 

Proceeding  in  an  admittedly  nonrigorous  fashion,  Ba.  (b4)  vas  treated  as 
a  qusidratic  expression  for  (valid,  hopefully,  at  and 

Subsequent  evaluation  of  this  expression  indicated  that  one  root 
gives  the  origin  value  correctly  and  the  other  root  gives  the  sonic  value 
correctly. 


Carrying  out  the  outlined  scheme  entailed  a  considerable  amoxait  of  work  in 
differentiating  and  expanding  the  expression.  Owing  to  its  length  (as 
well  as  its  limited  applicability)  the  resultant  eqmtion  is  not  Included 
here)  however,  a  plot  of  the  variation  of  the  two  roots  and  ^  is 
given  for  a  typicsG.  case.  Pig.  B2. 


The  proposed  scheme  was  not  euh>pted  in  view  of  the  following  considerations: 

a.  The  'value  did  not  approximate  with  sufficient  accuracy 

in  a  large  enough  neighboihood  of  the  sonic  point. 

b.  The  "switching  transi«its"  were  more  severe  in  some  causes  than 
erratic  functional  b^avior  due  to  stepping  too  close  to  o  . 


c.  Ihe  1.  value  was  sufficiently  sensitive  to  the  value  of  eT  to 

require  almost  the  same  degree  of  ^  refinement.  This  sensiti¬ 
vity  is  indicated  in  terms  of  emd  in  Jigs .  B3 

and  b4,  respectively.  Hie  practical  advantage  of  utilizing  the 
Z  expression  is  th\is  lost. 

d.  The  sheer  bulk  of  the  Z  expression  contributed  significantly 
to  computing  time  expended. 
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